Charge Fluctuations in the Single Electron Box 



O 

o 
o 

(N 

Oh' 

oo 






cd 



13 

o 
o 



> 

^t 

o 
o 
o 

-I— > 



I 

c 
o 
o 



x 



Georg Goppert and Hermann Grabert 

Fakultat fur Physik, Albert-Ludwigs-Universitat, 

Hermann-Herder-Strafle 3, D-79104 Freiburg, Germany 

(February 1, 2008) 

Quantum fluctuations of the charge in the single electron box are investigated. Based on a 
diagrammatic expansion we calculate the average island charge number and the effective charging 
energy in third order in the tunneling conductance. Near the degeneracy point where the energy 
of two charge states coincides, the perturbative approach fails, and we explicitly resum the leading 
logarithmic divergencies to all orders. The predictions for zero temperature are compared with 
Monte Carlo data and with recent renormalization group results. While good agreement between 
the third order result and numerical data justifies the perturbative approach in most of the parameter 
regime relevant experimentally, near the degeneracy point and at zero temperature the resummation 
is shown to be insufficient to describe strong tunneling effects quantitatively. We also determine 
the charge noise spectrum employing a projection operator technique. Former perturbative and 
semiclassical results are extended by the approach. 
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I. INTRODUCTION 

Tunneling of electrons in metallic nanostructures is 
strongly affected by Coulomb repulsion. Provided the 
screening length in the metallic films is small compared to 
tunneling barrier thickness and sample size, the Coulomb 
energy can be written in terms of a geometrical ca- 
pacitance C. The relevant energy scale of the system 
is the corresponding charging energy Ec = e 2 /2C |I]] 
that is the energy needed to charge the capacitance C 
by one electron. For temperatures well below this en- 
ergy, ksT <C Ec, tunneling onto a metallic island is 
exponentially suppressed. For weak tunneling strength, 
Gt <C Gk-, where Gt is the phenomenological tunneling 
conductance and Gk = e 2 /h the conductance quantum, 
systems are well described by the perturbative approach 
in a = Gt/Gk [PL[|- When the tunneling conductance 
becomes larger, higher orders in the perturbative series 
in a such as cotunneling ||] have to be included. Even at 
zero temperature these processes are not forbidden en- 
ergetically, therefore, higher order corrections are most 
pronounced at low temperatures where first order pro- 
cesses are exponentially suppressed. Due to the large 
number of terms in the perturbative series, the approach 
remains restricted to the first few orders in a and one is 
interested in the range of validity one could expect. At 
higher temperatures, fc^T 3> Ec, thermal fluctuations 
dominate and perturbation theory (PT) becomes more 
accurate. Therefore, to give a lower bound of the valid- 
ity of PT it is sufficient to consider the zero temperature 
case where PT is worst. In fact, PT at zero tempera- 
ture even diverges at the degeneracy point of the single 
electron box (SEB) showing that the range of validity 
depends on both, tunneling strength and applied gate 
voltage. 

While partial resummation techniques pHq] lead to 



a nonperturbative finite result even at the degeneracy 
point, they need an arbitrary cutoff that limits their use 
for direct comparison with experiments. Recently, renor- 
malization group (RG) ideas M have been used to re- 
move the cutoff yielding results that depend on param- 
eters measurable experimentally. On the other hand, a 
complete resummation of PT can be achieved in phase 
representation leading to a formally exact path integral 
formulation of single electron devices p0 , O |. Here, the 
phase is the conjugate operator to the island charge num- 
ber and can be related to the physical voltage drop across 
the junction. The functional can serve as starting point 
for analytical predictions in the semiclassical limit |lj-[l|] 
covering the range of high temperatures and/or large con- 
ductance. It also is the basis of numerical calculations 

In this paper we study the single electron box by sys- 
tematic diagrammatic techniques. So far PT for any 
temperature has been calculated to the first and the sec- 
ond orders in a ]24|| . Here, we determine the third order 
corrections at zero temperature and compare them with 
Monte-Carlo (MC) results §0|-§|,|| and recent RG data 
p6[ . Further, we explicitly address the vicinity of the de- 
generacy point where PT fails and derive a nonperturba- 
tive result by resummation of graphs contributing to the 
leading logarithmical divergencies. We also discuss the 
charge noise spectrum nonperturbatively by means of a 
projection operator technique. In general, the spectrum 
depends on the electronic bandwidth, however, for fre- 
quencies relevant experimentally only the high frequency 
cutoff characterizing the resolution of the measuring de- 
vice matters. 

The paper is organized as follows: In Sec. || we in- 
troduce the system Hamiltonian and the average island 
charge number. In Sec. [II we briefly recapitulate essen- 
tial results and diagrammatic rules of the perturbative 



expansion J2J]. The necessary changes and simplifica- 
tions in the zero temperature limit are given in Sec. IV 
where we exemplarily evaluate one graph of third order. 
In Sec. |V] we present the analytic result for the average 
charge number. The two state approximation and the 
resummation of the leading logarithmical divergencies at 



the degeneracy point are discussed in Sec. VI. In Sec. VII 



we compare the analytical findings for the average island 
charge number and the effective charging energy with MC 



data and recent RG results. Finally, in Sec. VIII we de 



termine the noise spectrum of the island charge number 



and conclude in Sec. IX 



II. SYSTEM AND MODEL HAMILTONIAN 

We consider a SEB consisting of a metallic grain that 
couples to a lead electrode via an oxide layer. The sep- 
aration permits tunneling of single electrons with the 
corresponding phcnomcnological tunneling conductance 
Gt- The geometrical capacitance between the grain an 
the lead reads Ct- Furthermore, a gate electrode is ca- 
pacitively coupled to the grain with gate capacitance C g . 
This setup is shown schematically in Fig. [fl, where the cir- 
cuit is biased by a gate voltage U g shifting the Coulomb 
energy continuously. We describe the SEB by the Hamil- 
tonian [51 



where 



H = Hq + Ht 



H — He + H„ 



represents the system in absence of tunneling. Here, 



H w ~ /2 e k<? a la a k<? 



/ j t qa a qo a qa 



(1) 



(2) 



(3) 



describes free Fermions where the quasiparticle creation 
and annihilation operators for transversal and spin quan- 
tum number a and longitudinal quantum number p = k,q 
for the island and the lead electrode, respectively, are 
denoted by a pa and a pa . Further, e pa are quasiparticle 
energies for the corresponding quantum numbers. 



H c = E c {n 



n g ? 



(4) 



is the Coulomb energy for n excess charges on the is- 
land biased by the dimensionless external voltage n g — 
CgUg/e. The charging energy 



H T = ^2 ( W a L*V A + H - 



kqa 



(6) 



with tkqcr the transition amplitude between states with 
quantum numbers ka and qa. The charge shift operator 
A accounts for the Coulomb energy and is related to the 
charge number operator n by the relation 

A^nA = n - 



1. 



(7) 
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FIG. 1. Circuit diagram of a single electron box. 

The excess charge number n can be expressed as a 
derivative of the system Hamiltonian with respect to n g . 
Accordingly, we find for the average island charge num- 
ber 

/ x 1 dlnZ 

( n ) = n 9 + 7^T -5 . (8) 



where 



2(3E C dn g 



Z = tr{exp(-/?#)} 



(9) 



is the partition function of the system. At low temper- 
atures and in the limit of vanishing tunneling conduc- 
tance the logarithm of the partition function reduces to 
the minimum of the electrostatic energy Ec(no — n g ) 2 
where no is the integer closest to n g . Hence, as a func- 
tion of the applied voltage U g the island charge number 
displays the well known Coulomb staircase (n) = no ob- 
served experimentally pTf] . Due to occupation of higher 
energy levels at finite temperatures the step function is 
smeared. Similarly, the Coulomb staircase is rounded 
by virtual occupation of higher charge levels caused by 
the finite tunneling conductance. Here, we restrict our- 
selves to zero temperature and discuss the influence of 
higher order tunneling processes on charge fluctuations. 
Because of the periodicity and symmetry of the partition 
function Z with respect to n g , it is sufficient to consider 
0<n g <i. 



En = 



2C 



(5) 



depends solely on the island capacitance C = Ct + C g . 
Spin and transversal quantum numbers are conserved 
during the tunneling process described by the tunneling 
Hamiltonian 



III. PERTURBATION EXPANSION AND 
DIAGRAMMATIC RULES 

In this section we briefly summarize the method in 
Ref. p4fl and give the diagrammatic rules used in the 
remainder. 



A. Perturbation Expansion 

Since the n g dependence of the partition function (0) 
arises from the charging energy only, we may put 



tre 



-8H 



tr qp e-P H o 



(10) 



Factorizing the exponential exp(—(3H) into a part 
exp(— /3-ffo) hi the absence of tunneling and an interac- 
tion part written as a series in the tunneling Hamiltonian 
Ht we get 

X fP P<y-m f&2 

Z= V](-l) m / da m / da m -i--- / da% 
m=o Jo Jn Jn 

tr{e-P H °H T (a m )...H T {a 1 )} 



tr qp e~^ 



(11) 



with the tunneling Hamiltonian in imaginary time inter- 
action picture 



rflo ir„„-Tffn 



H t (t) = e T "°H T e 



(12) 



Separating the trace in a Coulomb and a quasiparticle 
trace 



tr e- pHo . . . = tr c e _/3H ° tr aD e _/9H * 



iqp ' 



(13) 



one is left with multi-point correlation functions of tun- 
neling Hamiltonians in imaginary time 



Z= ^(-l) m / da m / dcc m -i... \ da! 
m _ n Jo Jo Jo 



tr c e' pH ° (H T (a m ) ■ ■ ■ H T (ai))o . 



where 



(X)o = 



tr qp {e-^X} 
tr qp e _/3ff< jp 



(14) 



(15) 



is the thermal quasiparticle average. Due to the Coulomb 
interaction these correlators do not decompose into a 
product of two-point correlators. However, inserting the 
explicit form (g) of Ht, the charge shift operators in 
interaction picture A(r) = exp(riJc , )Aexp(— tHq) com- 
mute with the quasiparticle operators and therefore may 
be factored out of the quasiparticle trace leading to 

Z=VV-l) m / da m dam-!... da\ 

rn=0 J ° J ° J ° 

E ••■ E im ci---c m 

fel <?1 &\ Cl fem 9m cr m Cm 

Y^ e - 2:^(^-^-0^ (16) 

( a lb<r m ( a ™)%i^ m ( a ™) ■ • ■ 4U ( a i) a ^i ( a i)) ■ 



Here, the Coulomb trace is explicitly represented as a 
sum over charge states labeled by n and 



E n = Ecin-ng)" 2 



(17) 



is the corresponding Coulomb energy. The charge shift 
operators in interaction picture lead to exponentials 
exp[aj(E n . +1 — E nj )] where the integers rij are defined 
by 



m = n, 



3-1 

i] = E & 



(18) 



with C,k — i labeling excess charge number increasing or 
decreasing processes. Due to charge conservation the £ 
sums in Eq. ( |16|) are constrained by 



2m 

3=1 



(19) 



Further, we have introduced the shorthand notation 
a + = a' and a - = a, respectively, and i = tk qa is a real 
averaged transmission amplitude. Now time difference 
variables (3j, (j = 1, ...,2m) between two subsequent 
operators are introduced and the cyclic invariance of the 
trace can be used to obtain 






m=0 



dpi 



j-oo I 2m \ 

J d(3 2m 5 IE&-0 



-V 2m /3,-E„. 



E E< 

Cl.---.C2m "=-oo 

(fic^t^t^^fE^ 



E ■•• E 

k 1 q 1 a 1 k 2m q2m.<?2r, 



(20) 



vi=l 



J=l 



Since the free quasiparticle Hamiltonian _ff qp is quadratic 
in the fermionic degrees of freedom, the expectations of 
quasiparticle operators ol a (T) and a pa (T) obey a Wick 
theorem and the average in (120) decomposes into a sum 
over pair products of two-point correlators 

e Cl(r\—T 2 )e plal 
(apU( T l)S2<X2( T 2))o = fri.-cAi.pAi,^ 1 + eCl /3e P1J1 ' 

(21) 

So far the result is valid for arbitrary numbers of tunnel- 
ing channels 



AA=£l. 



(22) 



In metallic junctions where the junction area is typi- 
cally much larger than the Fermi wavelength squared, 
J\f is very large. Experimentally the value is of the order 
J\f w 10 4 justifying limiting considerations. 




FIG. 2. Representative circle diagrams of a) first, b) sec- 
ond, and c) third order in the perturbative series. 



g / de k 



eke' 



Wk\/D 



1 



-Pe k 



(27) 



for k — 1, . . . , m, stemming from electron- hole pair Green 
functions (J23J). Since with respect to the auxiliary vari- 
able E the whole integrand is a product of energy de- 
nominators, the integral can be performed explicitly by 
means of contour integration. 



The 1/OV corrections for the SEB are considered explic- 
itly in Ref. P8fl confirming the validity of the approxima- 
tion. In leading order only the combination 

Y(t 2 - Ti ) 

= t 2 E («L.(^)a^(n))o(«-i(^)4 lCT (n))o 



fcl<Jlfc 2 <72CT 



g / de 



-H/£> 



1 - e -f>* 



,-(T2-Tl)t 



(23) 



of two two-point correlators contributes. Here, we have 
replaced the sums over k and q by energy integrals and 
have already performed one of them. We introduced the 
notation g = t 2 M pp' = a/An 2 where p and p' are densi- 
ties of states at the Fermi level for the island and the lead 
electrode, respectively. The function Y(t) is an electron- 
hole pair Green function where electron and hole are cre- 
ated in different electrodes and D is the electronic band- 
width. Representing the (5-function in Eq. ( |2C| ) in terms 
of an energy integral over an auxiliary variable E 



S[t) = tJ 



dEe 



-irE 



(24) 



one can perform all imaginary time integrals /3j, (j = 
1 . . . 1m) gaining energy denominators that are linear 
combinations of the auxiliary variable E as well as E nj , 
(j = 1 . . . 2m) and £&, (k = 1 . . . m). The coefficients of 
the linear combinations depend on the Wick decomposi- 
tion and the £-sums in (|20f). The remaining summations 
over pairs and Q's can be represented graphically by di- 
agrams, in terms of which the partition function reads 



1 r°° 

2W-00 ,.^ 

diagrams 



D(E) 



(25) 



In Fig. H representative circle diagrams of a) first, b) sec- 
ond, and c) third order in PT are shown. A diagram 
D(E) includes the charge sum X)^L-oo an< ^ eacn circle 
segment correspond to an energy denominator 



E nj + J2T=i Qjktk + iE 



(26) 



where 



jk 



1 for a time interval f3j in between of 
two vertices connected by a straight tunnelon line, and 
Ojfc = otherwise. The tunnelon lines represent energy 
integrals 



(1) (2) 

FIG. 3. Diagrams of first order in the perturbative series. 



One gains a sum over residua corresponding to poles at a 
certain circle segment. Since residua of higher order poles 
correspond to derivatives of the product of finite energy 
denominators with respect to E, we get graphs decorated 
by slashes. The decorations are placed on the finite seg- 
ments indicating a derivative of the corresponding en- 
ergy denominator with respect to the auxiliary variable 
E evaluated at the pole position considered. The simple 
form of the energy denominators allows us to perform 
the derivatives explicitly. We gain a higher power of 
the energy denominator l/£ q+1 times (— i) q ql for a q-iold 
derivative of l/£. This way the sum over residua can be 
represented by circle diagrams, where the divergent circle 
segments are omitted and the remaining finite segments 
are decorated with slashes indicating derivatives of the 
energy denominator evaluated at the pole position. Con- 
sidering a circle diagram with a pole of order r, one finds 
that when the pole segments are omitted, the graph de- 
composes into r pieces denoted by T q with q = 1, . . . , r. 
The partition function may then be arranged according 
to the order of poles 



r— 1 diag n— — oo 



-pE n 



E (V 

8=0 V 



/3 r 



1 



1)!/ 



.9=1 



w 



(28) 



where the sum over diagrams is restricted to those with 
poles of order r and the symbol [. . .]( s ) stands for the 
sum over all decorations (derivatives) between the brack- 
ets with s slashes. Further, the factor / is the number of 
identical subgroups in T q , q = 1, . . . , r. Introducing the 
quantity 



^ = E 

diag 



.9=1 



(») 



(29) 



we may write the partition function in the form 




n— — oo 



(s) 

Using analytical properties of Ur one may show that 



E 

s=0 



u p+s 



P + S 



(g^)' 



(31) 



This way one gets an effective Coulomb representation of 
the partition function 



*= E 



-/3(B„+A„) 



where 



oo .. 



r=l diag 



,3=1 



(r-1) 



(32) 



(33) 



is the energy correction to the Coulomb energy of state 
n that may be written 



A — V^ A( m ) 



(34) 



771=1 



,(m) 



where A„ is the contribution of order g r 



(m) 



B. Diagrammatic Rules for A„ 

In this subsection we give the diagrammatic rules for the 
energy corrections (34). The term of order g m is given 

by An and is composed of graphs containing a vertical 
line with m semicircles attached. Each semicircle cor- 
responds to a tunneling event and represents an energy 
integral 



9 



dc 



ce 



N/-D 



1 



-0e 



(35) 



stemming from the electron-hole pair Green function 
(|23|). Further each vertical line element contributes an 
energy denominator — 1/£ where £ is the excitation en- 
ergy during the corresponding intermediate state, which 
is the sum of the Coulomb energy difference 



Sp -^Wi+p ^n 



(36) 



and of all electron-hole pair excitation energies tj present 
in the intermediate state that are represented by the arcs 
that would be intersected by a horizontal line. There are 
two types of semicircles: inflected to the right or left, 
whereby the Coulomb state n is increased (decreased) by 



an arc to the right (left). For example in first order in g 
there are just two processes depicted in Fig. |[ Whereas 
the graph (1) increases the excess charge number by one 
and therefore represents 



A (i,i) = 



-9 



de- 



1 



M/d l 



(37) 



the graph (2) lowers the charge number and it's contribu- 
tion An ' is obtained from the previous one by replacing 
Ci -+€-!■ 



g)00(3 



(1) 



(2) (3) 



(4) 



g(PQ) 



(5) (6) 



(7) (8) 



g)4)0(3 

(9) (10) (11) (12) 

FIG. 4. Diagrams of second order. 



In second order, cf. Fig. ^, there are two semicircles 
dividing the vertical line into three parts. Each of them 
represents an energy denominator at the corresponding 
charging energy. The graphs (1) to (8) in Fig. [|are given 
by all possibilities to attach two semicircles to a vertical 
line. Using the rules given above the graph (1) in Fig. || 
correspond to 



A (2,l) = 



dc\ 



do 



eie -\*i\/D e2e -\ez\/D 



-Pei 1 



-/3e 2 



(6+£i) 2 6+ei+£2 



(38) 



representing two tunneling processes with three interme- 
diate energy denominators. Since an arc to the left hand 
side represents a tunneling process that lowers the charge 
number on the island, the contribution of graph (2) dif- 
fers from that of graph (1) by the replacement £2 — > 0. 
The first eight graphs in Fig. || can be generated easily by 
the rules given. The graphs (9) to (12), however, differ 
by the prolongation of the interior arcs across the ver- 
tical line. These "insertions" represent separate graphs, 
i. e. each insertion represents a graph of lower order mul- 
tiplied to the main graph. However, the vertical line of 
the main graph has two peaces with the same energy 
above and below the insertion and therefore the energy 



denominator is squared. Moreover, each insertion carries 
a factor (—1). 



FIG. 5. Representants of distinct graphs of third order 
without insertions. 



Analogously, with two insertions the main graph has a cu- 
bic energy denominator, see for example the graph (16) 
in Fig. pi Using these rules, the graph (11) in Fig. Elleads 
to the energy correction 



At 2 - 11 ) 



dci 



do 



eie -Nd/^ e2e -|e 2 |/D 



— OO 

1 



-/3e x l 



-0e 2 



(£- 



ei) 2 6 



£2 



(39) 



that factorizes into (—1) times the graph (2) in Fig. y 
with the energy denominator squared, multiplied with 
the graph (1) in Fig. [|. These insertions are shorthand 
notations for graphs with decorations that stem from a 
pole of higher order in the energy denominator multi- 
plied with a lower order graph. Therefore the higher 
order denominators and the factors (—1) correspond to 
derivatives with respect to the auxiliary variable E. At 
zero temperature these graphs are related to terms in the 
Rayleigh-Schrodinger perturbative expansion stemming 
from normalization 1231 ■ 




(13) (14) (15) (16) 

FIG. 6. Representants of distinct graphs of third order with 
insertions. 



In the high temperature limit, or equivalently in the 
limit Ec — > 0, the Wick theorem implies that only "con- 
nected" graphs appear in the exponent. The "connected" 
graphs are those where all processes occur in one channel. 
In order g m , m > 1, these "connected" graphs, however, 
give only 1/M corrections that we have omitted. Hence, 
in the high temperature limit only the first order graphs 
survive and all higher order terms have to cancel. For 
the second order terms it turns out that each column 
in Fig. H cancels. For example, the contribution of the 
first column, graphs (1), (5), and (9), leads at vanishing 
charging energy to 



A(2;l,5,9) 



E c =0 



r°o eie -M/D tze-M/D 
dt\ I de 2 - 



1 



1 - e-l 3 ^ 1 - e~P e 
1 1 



q(ei+e 2 ) ei(ei + e 2 )e 2 qe 2 



0. 



(40) 



Higher order diagrams cancel likewise, and therefore in 
the limit of a large channel number the perturbative ex- 
pansion can be truncated after the first term at high 
temperatures where the charging energy Ec is negligi- 
ble compared to the thermal energy fc^T. 

The application of the diagrammatic rules to higher 
order graphs is obvious and we just present the contribu- 
tion corresponding to graph (7) in Fig. |j| 



A (3,7) 



de\ 



do-, 



des 



[(fi + ei)(6 + ei + e 2 ) 



e ie -l'il/-P e 2 e-l £2 l/ D e 3 e-^ D 
1 - e-^i 1 - e-l 3 ^ 1 - e-P's 
x (6 + ei + e 2 + e 3 )(6 + ei + e 3 )(£i + e 3 )] _1 (41) 

as an example of a third order term. Here, three tunnel- 
ing processes occur and we have to deal with five energy 
denominators. In this graph the excess charge number 
is raised three times and then lowered stepwise. The 
energy denominators include the quasiparticle excitation 
energies of arcs that would be intersected by a horizon- 
tal line. While the creation and annihilation times aj, 
(j = 1, . . . , 2m) of these quasiparticle excitations are al- 
ready integrated out, the order of creation and annihila- 
tion of different excitations is crucial. All types of dia- 
grams of third order PT without insertions are depicted 
in Fig. H and the graphs with one or two insertions are 
shown in Fig. pi Here, we have omitted different combi- 
nations of semicircles to the left and right but displayed 
only one representant. Hence, each graph in Figs. and 
o stands for 8 different left-right combinations of the 3 



semicircles. Further, the diagrams (7) and (8) in Fig. |] 
and (11) and (12) in Fig. g are topologically distinct from 
graphs reflected at a horizontal line, but obviously, they 
lead to identical contributions and we have to count them 
twice. 

In summary, we have 80 topologically different graphs 
without insertions, 64 with one and 16 graphs with two 
insertions leading to 160 topologically different graphs of 
third order. All contributions of third order are readily 
evaluated using the rules given above. In general, there 
are always 2 m left-right combinations for one represen- 
tant of order m. Moreover, the number of representants 
exceeds the (2m — 1)!! possibilities to arrange m semi- 
circles along the vertical line, simply as a consequence 
of the summation over insertions. Hence, the number 
of graphs of order m exceeds (2m)!/m!, i.e. grows faster 
than the factorial. This rough estimation leads to more 
than 120 graphs in third and more than 1680 graphs in 
fourth order, showing that a reasonable treatment is lim- 
ited to third order. To obtain higher order or nonper- 
turbative results, one has to limit oneself to partial sum- 
mations of graphs including the essential contributions. 
Unfortunately, in the infinite cutoff limit D — > oo, each 
graph of the perturbative series represents a diverging in- 
tegral and only the full sum in each order remains finite 
p4|. Hence, partial summations of higher order graphs 
need an artificial cutoff complicating a direct comparison 
with experimental findings. Consequently, within PT, 
the systematical treatment of higher orders is the only 
tool to get results directly comparable with experiments. 
To proceed we discuss two general simplifications, valid 
for all orders. 



C. Reflected Graphs 

The analytical form of the integrals leads to general con- 
sequences for energy corrections. First we note that the 
particular form of the charging energy (|j) implies 



E n {n g ) = E n+ i(n g + 1) = E_ n (-n g ) 



(42) 



and the corresponding energy differences in the denomi- 
nators read 

C± P = E n±p (n g ) - E„(n g ) = E c [p 2 ± 2p(n - n g )] (43) 

which depends, on the difference (n — n g ) only. 





JP ^ J -P 

FIG. 7. Symmetry of two vertically reflected graphs. 



Hence, the corrections of order m may be written in the 
form 



AlrHn g )=g m Ecf m (n~n g ). 



(44) 



Since contributions of order m include always m integrals 
with 2m — 1 energy denominators, we gain by measuring 
all energies in units of Ec a single factor Eq , and in the 
limit of an infinite bandwidth, D — > oo, the functions 
fm depend on flEc only. Further, a reflection of a given 
graph on the vertical axis leads to the same contribution 
with the replacement E n ± p — > E n ^ P) cf. Fig. ffl. Equiva- 
lently, by virtue of Eq. (H3), one can replace (n — n g ) by 
— (n — rig). Since the sum over diagrams includes all left- 
right permutations of arcs including pairs of vertically 
reflected graphs, one may write 



fm(u) = g m (u) + g m (-u), 



(45) 



where the contribution of g m solely includes topologically 
different graphs where one arc is held fixed. 



D. Insertions 

A further simplification arises for graphs with insertions. 
It turns out that they factorize into a host graph, with 
energy denominator squared, and an insertion contribu- 
tion. Since we have to sum over all possible left-right 
configurations, the full lower order contribution may be 
inserted. In Fig. |SJ a) we replace a first order insertion 
and its vertically reflected companion by a circle repre- 
senting a multiplication with the full first order contribu- 
tion A« . Likewise, insertions of order k and all possible 
left-right combinations lead to a multiplication with the 
full order k contribution A n , schematically depicted for 
k = 2 by the square in Fig. b). The double arrow 
in this figure represents the sum over all possible left- 
right arrangements of semicircles belonging to the inser- 
tion. Therefore, all graphs belonging to the representants 
(13), (14), and (15) in Fig. | lead to 
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fo 



(46) 



with u = n — n g . This is a multiplication of three fac- 
tors, A„ ' and the host graph (1) in Fig. g with the energy 
denominator squared, and (—1). Additionally, we have 
added the contribution of the vertically reflected graphs 
which is of the same form with the replacement u — > — u. 
Likewise, all graphs belonging to the representant (16) in 
Fig. H lead to 



A( 3 < 16 > = -g*Ey x {u? 
+ {u — > —u) , 



de 



-WD 



1 - e-P* (& + e) 3 



(47) 



which is a multiplication of I A„ j and the host graph 

(1) in Fig. with the energy denominator cubed. This 
procedure holds for insertions of all order and is used in 
Sec. |v| to resum the leading logarithmic divergencies in 
the two-state approximation. 



a) 

b) 
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FIG. 8. a) full first and b) full second order insertions. 



IV. ZERO TEMPERATURE LIMIT AND THIRD 
ORDER 

In this section we perform the zero temperature limit and 
present results of third order. 



A. Zero Temperature Limit 

Generally, at zero temperature the partition function of 
a system with nondegenerate discrete energy levels de- 
pends only on the ground state energy 



E 



-/3(£„+A„) 



-pe 



(48) 



Due to the restriction of the gate voltage to the range 
< n g <i, the ground state charge number is n = and 
the ground state energy reads £ = E + A . Therefore, 
the energy differences in the denominators read 



£±p = E± p -E = E c (p 2 T 2pn g ) 



(49) 



and the functions f m and g m depend on n g only. The 
average island charge number reduces to 



1 d£ 

(n) = n„ — — 

W 9 2E C dn a 



(50) 



Bose factors in the integrands restrict the integration ( J35| ) 
to positive energies 
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de- 



,-\A/d 



1 



= -/3e 



fJ 



deee~*' D 



(51) 



which facilitates the integration, because there are no 
poles from the energy denominator (pq) that need to be 
taken care of, meaning there are no real excitations. 



B. First and Second Order Results 

The first and second order calculations were already pre- 
sented in Ref. 24 1. The first order with an exponential 
cutoff D leads to 



9i(u) = -(1 + 2m) ln(l + 2u) 



D/E c + ln(D/E c ) + 7, 

(52) 



where 7 = 0.577. . . is Euler's constant. In the infinite 
cutoff limit this term diverges, but the derivative with 
respect to n g remains finite. When we restrict ourselves 
to first order in g, we may omit this divergence, however, 
in combination with higher order terms in the perturba- 
tive series, i.e. as lower order insertion, we have to use 
the full expression (|52|) ■ Whereas the first order contribu- 
tion diverges, the full sum of graphs of each higher order 
m > 1 remains finite (24j . The second order contribution 
was also calculated previously leading to 
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4(1 -u) In 



/ l-2u 

U(i-«) 



-(5-8M + 4w 2 )Li 2 



3-2u 



A(l-u)J ' 
where Li 2 (w) is the dilogarithm function J29| 



(53) 



C. Third Order 

Here, we motivate the calculation of the third order con- 
tribution exemplarily for the graph (7) in Fig. |5j. The full 
analytic expression 33(1*) is presented in the Appendix. 

Since the integral over the entire sum of 160 integrands 
remains finite, we may omit the cutoff. However, to cal- 
culate the integral analytically we have to separate the 
whole expression into tractable parts. Each integral will 
be divergent and we introduce a sharp high energy cutoff 
D. After integration we expand the expressions with re- 
spect to D/Eq — ► 00. There are divergent terms in each 
expression but the sum of divergencies of all graphs has 
to vanish. This cancellation serves as a useful, nontriv- 
ial test of our calculation. We have to deal with eight 
different types of integrands without insertions depicted 
in Fig. and six different types with insertions in Fig. o 
(the graphs (13) to (15) are merged to a single graph 
with the whole second order contribution inserted). We 
exemplarily proceed with graph (7) in Fig. ra leading to 
an integral of type 



D pD pD 

de x I de 2 / de 3 e 1 e 2 e 3 



[(«1 + Ci)(«2 + e a + e 2 )(K3 + £l + £2 + £3) 
x(K 4 + ei+e3)(K5 + e 3 )] _1 . (54) 

where the Kj stand for excitation energies of the form 
(p6|). The full contribution of this representant consists 
of all possible left-right arrangements of the semicircles 
leading to 



A 



(3,7) 



-9 



} j hrJKl, K2, K3, K4., K5), 



(55) 



where the sum runs over all allowed combinations 
of energy differences. These are h-j (^1,^2,^3^2,^1), 
hr(£ifl&&,ti), hrfaAZ-iM-i), h 7 (£ 1 &&,0,Z-i), 



and all terms with n 



The integrals may be per- 



formed by splitting denominators into partial fractions. 
Using 

1 1 1 



(a + e 2 ){b + e 2 ) 



{a-b)(b + e 2 ) (a-b){a + e 2 ) 

(56) 



where a = k 2 + e\ and b = K3 + £1 + £3, we are able to 
perform the e 2 integral leading to a logarithm function 



a ln(a + D) - a ln(o) - D. (57) 



de 2 — ; 

a + e 2 



The new denominator l/(a — 6) on the rhs of Eq. Q56 
has an artificial pole at £3 = n 2 — K3. Since the integrand 
is analytic in the integration region, the sum of all pole 
contributions has to cancel in the threefold integral d54). 
However, each pole contribution depends on the contour 
of integration and we have to specify the contour and 
use the same for all integrals. Next we consider the €1 
integration. In the numerator now appears the logarithm 
function from the previous integration, ln(«;2 + £i), that 
diverges for t\ — > when k 2 = and we temporarily in- 
troduce a lower integration limit. With a decomposition 
of the form (p6|), where £2 is replaced by £1 and the con- 
stants read a = K\ and b = K4 + £3, we split the remaining 
fractions in (p4| ) and perform the integral in terms of the 
dilogarithm function [E9J 



Li2(.s) = dz 

Jo 



,ln(l-z') 



(58) 



The third integration can then be performed using the 
trilogarithm function Li 3 (z), where the general polylog- 
arithm functions are defined by |29] 



Li n (z) — I dz 



,Li n _i(z') 



(59) 



All terms emerging can be expressed by trilogarithms, 
dilogarithms, logarithms, and rational functions of n g . 



Here, the arguments of the transcendental functions are 
rational expressions of energy differences. Since each in- 
tegration increases the "order" of transcendental func- 
tions at most by one, the transcendental terms are of the 
form Life In obeying k + l < 3. Here In stands for a prod- 
uct of I logarithms of possibly different arguments. The 
"order" thereby characterizes the transcendental func- 
tion for large arguments: e.g. Life ~ In is of order k. 
We find in order m of the perturbative series transcen- 
dental terms of the form (Li m ) fcm • • • (Li2) fc2 In :i where 
Y^JLi jkj < rn. In principle, the integrals in PT lead 
in all orders to analytically known functions but there 
are practical restrictions, in particular, since the inte- 
grals cannot be evaluated straightforwardly by tools like 
Mathematica because one has to take care of the integra- 
tion contours and pole contributions explicitly. In view 
of the length of the analytical result we present 53 (u) in 
the Appendix. 

The result is in terms of complicated analytical func- 
tions that have to be calculated numerically. Therefore, 
one could think of a direct numerical evaluation of the 
three fold integrals, but there are serious numerical prob- 
lems. First, only the full sum of the 160 integrals is finite 
and one has to use a huge integrand. Moreover, the inte- 
grand is not symmetric in the three integration variables, 
in particular, there are integration directions where the 
integrand leads to diverging positive and negative con- 
tributions. For a numerical study one has to symmetrize 
the integrand which enlarges it by a factor of 6. Second, 
in spite of the smoothness and analyticity of the inte- 
grand and the absence of poles in the integration region, 
the integrand contains oscillatory parts so that standard 
numerical routines fail. We used a statistical integration 
routine to check the analytical predictions for selected 
parameter values. Typically, the numerical evaluation of 
a single point with a few per cent accuracy took over a 
week of CPU on a SGI Origin 200. Therefore, a numer- 
ical evaluation of the integral in third order is hopeless, 
in particular for delicate high precision studies regarding 
the limiting behavior near the degeneracy point or calcu- 
lations of second order derivatives needed to determine 
the charging energy. 



V. AVERAGE CHARGE NUMBER 

The analytical result for the average island charge num- 
ber (n) at zero temperature and in first order in g was 
calculated previously |5pQ] 



, l + 2n g 

Wi -.9 In- —^ 

1 — 2n„ 



(60) 



and can be readily obtained by using Eq. (15(1) with the 
energy correction (p2). The function is well behaved ex- 
cept at the degeneracy point n g =| where it exhibits a 



logarithmic divergence (n) 1 ~ —glnS, with 8 = \ — n g . 
Therefore, the range of validity of the perturbative series 
at zero temperature is strictly restricted to n g <a. The 
contribution of second order in g presented in Ref. |3l| ] 
reads 
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4tH 2 ( \-2n g 
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(61) 



16(1 + 2n„ - 2n 2 ) 
+ ,„ - w ! . n 9 { ln(l - 2n 9 ) 



(3-2n g )(l + 2n g ) 
+2(1 -%) 



ln 2^ 1-2% 



4(1-%) 



2Li 2 ' 3 " 2 ^ 



8(1~%) 



(l-2n 9 )(3-2n 9 ) 



ln(4(l - n g )) 



4(1-%) 



- (% -> -%) 



Here, (n fl 



— n„) stands for the same sum of terms 



with n„ replaced by 



showing explicitly the asym- 



metry of (n) with respect to the applied voltage U g . At 
the degeneracy point, a leading logarithmic divergency 
of (11)2 ~ — 2g 2 \r? 5 appears. This indicates that near 
the degeneracy point g\n8 is the effective expansion pa- 
rameter so that the larger g the smaller is the range of 
n g where finite order PT suffices. The analytical result 
of third order is not given explicitly here but can easily 
be calculated by differentiating the expression A with 
respect to n g , cf. (p0|). Also the third order term shows 
a logarithmic divergency at the degeneracy point lead- 
ing together with the lower order contributions to the 
asymptotic expansion 



(n) = ag 2 + bg 3 - (g + 6g 2 + eg 3 ) In 5 

-{2g 2 + 24 5 3 ) In 2 S - Ag 3 In 3 6 + O (6, g 4 

where the coefficients a and c read 



(62) 



VI. DEGENERACY POINT 




FIG. 9. Full ground state energy diagram 
non-crossing approximation of the two-state model. 



the 



A further analysis of the perturbative series in the limit 
of n g — >§ shows that the leading logarithmic divergen- 
cies ( pq ) stem from diagrams including only charge states 
11 = and 1 . This is a direct consequence of the degen- 
eracy of these charge states at n g =\. The two-state ap- 
proximation limits the perturbative series to graphs that 
contain the charge states n — 0, 1 only. This model was 
shown to be equivalent to an anisotropic multi channel 
Kondo Hamiltonian where EqS corresponds to the mag- 
netic field and g to the exchange integral || . Considering 
the leading order logarithmic divergencies we find that 
crossed diagrams, like graph (6) in Fig. [^ with the mid- 
dle semicircle to the left, do not contribute. Therefore, 
we may restrict ourselves to noncrossing diagrams. It is 
then possible to write the ground state energy as 



An 



deeiq(e) 



(66) 



graphically represented in Fig. H. The generalized energy 
denominator iq(e) corresponds to the bold line with in- 
dex 1 and is determined by a Dyson equation graphically 
represented in Fig. |lC|. For convenience, we have rotated 
the graphs by 90 degrees and an upper (lower) semicircle 
increases (lowers) the charge number. Further, since we 
sum a subset of graphs, we need to introduce a cutoff D, 
which for convenience is chosen as a sharp cutoff. 



a = 3--7r 2 -^ln(3) + 31n 2 (3) + 6Li 2 ( '- 



3 2 

= -9.7726..., (63) 

c = 21 + — 7r 2 + 91n(3)-61n 2 (3)-12Li 2 (| 
= 65.462 .... (64) 



+ 



+ 



In view of its length the coefficient b is given only nu- 
merically, b = —70.546..., but it can be readily obtained 
from the analytical result of A ' in the Appendix. The 
leading order logarithmic terms in Eq. (p2) read 



(n) g\n5-2g 2 ln 2 S 



3 ^ 3 S. 



4g 3 In 



(65) 



They are related to diagrams that contain only the charge 
states n = and 1. Therefore, before comparing our find- 
ings with numerical results, we consider a two-state ap- 
proximation restricted to these degenerate charge states. 



© = — + 



o 



FIG. 10. Dyson equation in the non-crossing approxima- 
tion of the two-state model. 
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In the Dyson equation depicted in Fig. [LO] the thin line 
with index corresponds to the bare energy denominator 
— 1/e of the charge state where e is the energy variable. 
Analogously, the thin line with index 1 correspond to 
the bare energy denominator — l/(£i + e) of charge state 
1. On the other hand, the bold line with index (1) 
represents the dressed n = (n = 1) propagator Fq(c) 
(-Fi(e)) arising from the insertion of interaction vertices 
according to the first (second) line in Fig. [l(| Hereby 
the interaction vertex represented by a circle with in- 
dex (1) is composed of a semicircle and an insertion 
according to the third (fourth) line in Fig. 10. These 



diagrams are evaluated by employing the rules given in 
Sec. |nj. Each part of an interaction vertex contains an 
integral g L de'e' but the two pieces differ in the explicit 
meaning of d: While the insertion is just a multiplica- 
tion where the propagator inside does not depend on the 
energy variable e of the legs of the vertex, the semicircle 
propagator takes into account all excitation energies and 
therefore depends on e + e'. Due to the restriction of the 
perturbative series to the charge states n = and n = 1 , 
the interaction vertices contain only two parts shown in 
Fig. fUJ and the Dyson equation generates all diagrams 
in the non-crossing approximation. However, to be con- 
sistent with the truncation to two charge states, one has 
to introduce a cutoff D that restricts electron-hole pair 
excitations to energies lower than the next charge state. 
With higher energy excitations being eliminated the pa- 
rameters must be interpreted as effective renormalized 
quantities. 



a) 
b) 

c) 
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FIG. 11. Graphs of a) first, b) second, and c) third order 
in the two-state model generated by the Dyson equation in 
Fig. 0. 



One way to find the renormalization of the parameters 
is a comparison of the limiting divergent behavior for 
n g ~~ *h with the full PT result. This comparison leads 
to a series expansion of the renormalized parameters g* 



and Eq in terms of the bare conductance g. 

Iterating the Dyson equation we generate graphs in- 
cluding only charge states n = 0, 1. In Fig. 11 we depict 



all graphs of a) first, b) second, and c) third order in g 
obtained this way. However, not all of these graphs con- 



tribute to the leading behavior near n g 



We find that 



only graphs of the form shown in Fig. [13 are responsi- 
ble for the leading logarithmic divergencies in order g m . 
The diagram a) shows m — 1 downward arcs leading to 
the asymptotic behavior 



A 



(m,0) 



-g m 6ln m d 



(67) 



for S — ► 0. The diagram b) shows m — k — 1 arcs lower- 
ing the charge state and one multi-insertion represented 
as an insertion with a slash with label k (k < m) which 
stands for the full result in the two-state approximation 
in order k. Assuming that the order k result has the 
leading order behavior 



;w 



\{2g\n8) k 5 



as is apparent for k = 1,2 and 3 form Eq. 
graphs asymptotically behave as 



(68) 
5|), these 
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6ln m S. (69) 



Interchange of the multi-insertion and the arcs leads 
to topologically different contributions that have to be 
counted separately. This leads to a factor (m — k) can- 
celing the denominator in Eq. (J69|). The sum over all /c's 
starting with k = as represented by graph a) in Fig. n3 
and from fc = ltofc = m— 1 corresponding to graph b) 
leads to an asymptotic behavior of the form 



A, 



(m) 



£4' 

fe=0 



(m,k) 



-1(2.9 In SY 



(70) 



which proves the assumption (|68|) for all k by induction. 
For the average island charge number (pCl) we obtain by 
summing over all orders in g 



-glaS 



l-2glnS 



(71) 



This result was originally obtained by Matveev g] with 
RG techniques. While this result does not depend explic- 
itly on the cutoff D, it contains renormalized parameters 
g and 5 depending on the tunneling strength. The leading 
logarithmic divergencies up to third order coincide with 
those given in Eq. (pa), however, nonleading divergencies 
from Eq. (p3) are missing. They can be included by using 
renormalized parameters g* and S* in Eq. (|7l|). Since the 
energy difference between the two charge states is Eq5, 
the renormalization of the charging energy may be writ- 
ten as EcS* . The renormalized quantities are then found 
to read 
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and 



g* =g[l + 6.g + c 2 g 2 + 0(g 3 )] 



5* = 5 [1 + d l9 + d 2 g 2 + (D(g 3 )] 



(72) 



(73) 



where the renormalization factors are given as series in 
powers of g with coefficients that may be expressed in 
terms of the expansion coefficients in Eqs. (p2|) (J64|) . We 
have c 2 = -4a - c = -26.372 . . ., d\ = -a = 9.7726 . . ., 
and c?2 = 6a + a 2 /2 — b = 59.662 .... As a consequence 
of the rapidly increasing series coefficients the result J7l| ) 
with (J7^) and ( f73| ) does not lead to meaningful results in 
the vicinity of n g =\ except for very small g. Therefore, 
one has to resum also nonleading divergencies in the two- 
state approximation to describe the behavior near n g =J. 
This has not been considered so far but is the aim of fu- 
ture work. 




FIG. 12. Diagrams of order g r 
logarithmic divergencies. 



m-k-1 

contributing to the leading 



VII. DISCUSSION OF RESULTS 

In this section we compare our analytical results with nu- 
merical data and estimate the range of validity of various 
orders of PT. While various QMC studies of the single 
electron box are available [P0|-p2 25 1, only the data in 
Ref. p5| determine the island charge number for finite 
gate voltages and very low temperatures so that they 
can be compared with zero temperature predictions. Fur- 
ther, we confront our results with the findings of a recent 
renormalization group study J9J]. The comparison of first, 
second, and third order PT with QMC data |25[ and RG 
results H for a = 2, 5, and 10 is shown in Fig. |13|. Not to 
overload the graph for a — 2 the RG results are omitted 
but they coincide with third order PT. We give results 
here in terms of the dimcnsionlcss tunneling conductance 
a = Gt/Gk = 47r 2 g. We find good agreement for gate 
voltages near zero, but for n g — >J the analytic result di- 
verges. As discussed above the range of validity of PT 
shrinks with increasing gate voltage. We find that third 



order PT in a remains valid with errors below 4% up to 
n g w 0.495 for dimensionless conductance a — 2, up to 
n g ~ 0.45 for a = 5, and up to n g rs 0.4 for a = 10. In 
these parameter intervals PT agrees both with QMC and 
RG data. Since for n g — 0.45 the charging energies for 
n = and n = 1 differ only by 0.1 Eq, deviations from 
the third order result in a can be observed only for tem- 
peratures well below Ec/lOks even at n g = §. Hence, at 
finite temperature the range of validity of PT increases. 
Further, Fig. |l3| shows that the resummation of the lead- 
does not suffice to describe 
, Subleading logarithms are 



Fig.^ 
ing logarithmic terms 
the behavior near n g 
important to obtain quantitatively meaningful results in 
the strong tunneling regime. We remark that the inclu- 
sion of subleading logarithmic terms of low orders only, in 
terms of the renormalized parameters g* and 5* defined 
in Eqs. (f72|) and (f73|), docs not improve the agreement, 
rather one has to consistently resum the nonleading log- 
arithmic terms to all orders. 
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FIG. 13. The average electron number (n) as a function of 
the dimensionless voltage n g is shown in first, second, and 
third order perturbation theory in a, and compared with 
QMC data Q and RG results |]. The result @ is also 
shown as a dotted line. 
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Since for n g — > the QMC and RG data perfectly 
coincide with the perturbative results, we used a more 
sensitive quantity to determine the range of validity of 
PT in this limit: For small external voltages, the average 
island charge grows linearly as 



(Q) = e(n) = C*U g 



(74) 



where C* is an effective capacitance of the box. In the 
absence of Coulomb blockade effects C* = C, while for 
strong Coulomb blockade, i.e., in the limit of vanishing 
tunneling conductance, C* = 0. It is thus natural to 
characterize the strength of the Coulomb blockade effect 
by an effective charging energy Eq defined by |2(J 



Eh 
E c 



1 



C* 

~c 



1 



d(n) 



dn a 



(75) 



Tl„=0 



One may view Eq as the effective rounding of the energy 
parabolas at n g — in the presence of tunneling. The 
perturbative series gives 



g = l-4 5 + V-5.9 3 + 0(3 4 ), 



(76) 



where 



^8^_32 + | ln(2) _ 161n2(2) _ 8Li2 g 
= 5.066... 



(77) 



and B = 1.457... is given here only numerically in view 
of the length of the analytical expression. 
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FIG. 14. Effective charging energy as a function of the di- 
mensionless conduction a — GtJGk- Perturbative findings 
are compared with QMC data |2C1-|22,|25[ and RG results M. 



In Fig. |lj we compare our predictions for the effective 
charging energy with QMC data by Wang, Egger, and 
Grabert @ (WEG), Hofstetter and Zwerger @ (HZ), 
Herrero, Schon, and Zaikin p^| (HSZ), and Goppert et al. 
H (GGPS) and RG results by Konig and Schoeller [§. 
We find good agreement with PT up to a ~ 8 for second 
order PT, while third order PT extends toa« 16. Dis- 
crepancies between these QMC studies only arise outside 
of the range of validity of third order PT. 



VIII. FLUCTUATIONS OF THE AVERAGE 
CHARGE NUMBER 

In this section we study the steady state fluctuations 
of the average charge number (n) . We are mainly inter- 
ested in the variance (Sn 2 ) with Sn = n — (n). Since 
(n — n g ) 2 = dHc/dEc, one may calculate this quan- 
tity in terms of a derivative of the free energy. At zero 
temperature, the free energy to first order PT is given 
by Ec\gi(n Q ) + gi(—n g )] with the function g\ defined in 
Eq. (|5^). However, this result would lead to a logarith- 
mic divergence of (Sn 2 ) for large electronic bandwidth 
D. The same divergence occurs also at finite tempera- 
tures, cf. Ref. ]24|] . To explore where this cutoff depen- 
dence comes from, we consider the correlation function 
(Sn(i)Sn(0)) . Using PT this quantity is well defined for 
finite t but diverges in the limit t — > 0. This divergence is 
found to result from al/w tail for large 10 in the fluctu- 
ation spectrum arising from the coupling to high energy 
electron-hole excitations. While this tail is cut off by the 
finite bandwidth D, in real experiments there is a cutoff 
from the time resolution of the measurement device at 
a certain frequency uj m well below D/h. To avoid these 
cutoff effects, we discuss only the noise spectrum 



S„(u>) 



dte iut C{t) 



(78) 



that is the Fourier transform of the symmetrized corre- 
lation function 

C{t) = -(Sn(t)Sn(0) + Sn(0)6n(t)) = {5n(t)5n(0)) s . 

(79) 

While S n (u>) also depends on the electronic bandwidth, 
for frequencies well below the electronic cutoff and not 
too low temperatures the spectrum becomes independent 
of the bandwidth. Here, 8n(t) = exp(iht)n — (n) is the 
fluctuation of the average charge number in the Heisen- 
berg representation, where the time evolution is governed 
by the Liouville operator LA" = j- [H, X] . 

Since a direct perturbative expansion of C(t) results 
effectively in a short time expansion, it is insufficient to 
determine the low frequency behavior of S n (u>) properly. 
Here, we use projection operator techniques |32j to de- 
rive a formally exact integral equation for the dynamics of 
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C{t) and then expand the kernel in powers of the dimen- 
sionless conductance g. The adequate projector reads 



px = wfW' 



(80) 



which fulfills the requirement P 2 = P. Applying the time 
evolution operator one finds the relation 



Pe lLt P = C(t)P , 



(81) 



with the normalized correlator C(t) = C(t)/(5n 2 ). Fur- 
ther, the time derivative of the correlator may be written 
in the form 



C(t)P = Pe lLt iLP. 



Using the operator identity 
ft 



iV L * = Pe* Lt P - 



dsPe^ s PiLPe 



iPLP(t-s) 



(82) 



(83) 



with P = 1 — P, Eq. ( JS2| ) may be written as 
C{t)P = Pe lLt PiLP 

dsPe ihs PiLPe IpLp{t - s \LP . 



(84) 



Inserting the identity ( |3l"| ) we get the exact evolution 
equation 



C(t) = ^C(t) 



ds<f>(t - s)C(s) 



(dn 2 ) 
with the memory kernel 

(nn r {t)) s (5nn) s (5nn r (t)) s 



4>(t) 



(5n 2 ) (Sn 2 ) (Sn 2 ) ' 

Here, the reduced time evolution of h is given by 
h r (t) = e ipLpt n. 



(85) 



(86) 



(87) 



The linear equation (85) can be solved by means of a 
Laplace transformation yielding 



C(z) = 



C(0) 



z — {Snn) s /{8n 2 ) + 4>{z) 
which is related to the noise spectrum by 



(88) 



S n (oj)- 



/OC /«OC 

dte iuJt C(t) = 2Re / die 1 
-oo JO 



l c{t) 



= 2Rc lim 



dte- zt C(t) = 2Rc lim C(z) . (89) 



To evaluate the right hand side of Eq. (Bq) we first deter- 
mine the initial value C(0) = (n 2 ) — (npin second order 
in the tunneling Hamiltonian by means of a straightfor- 
ward expansion of the density matrix using the imaginary 



time methods in Sec. Ill A. Since n and n 2 are diagonal 
in charge representation, both averages may be written 
in the form 



<X> = itr{e-^A} 



hi* 



-0H O 



dT'H T {T') 



f dr f dr" H t {t')H t {t")\xX (90) 

dT"Y(r' -T")e- {T '- T "^x) , 

I c 



-§<*+E 



C=± 



dr' 



with the free partition function Z = tr exp(— (3H ) and 
the Coulomb energy differences ^ that were introduced 
in Eq. (|36|). Here, we have decomposed the t race i nto a 
Coulomb and a quasiparticle trace as in Sec. Ill A, and 
have introduced the free Coulomb average 



(x )c = £ e-^x / y: 



-0E n 



(91) 



Inserting the representation (|23| ) for the electron-hole 
pair propagator into Eq. ( |9fj| ) , the time integrals are read- 
ily evaluated leading to an energy integral that may be 
solved by contour integration, cf. Ref. |2J. For the av- 
erage charge number squared we find 



(n 2 ) = (n 2 ) c (l + (^ 1 ) c ) 



(n 2 K 1 + 2nK 2 + K 3 ) c + 0(g 2 ). 



(92) 



The first term {K\)c stems from the expansion of the 
denominator Z where 

K x = 2irglm ]T [fi(iE c /v) - /i(i&/")] (93) 

C=± 

and v — 2ir/{3. A divergent part for D — > oo is omit- 
ted since it is canceled by a corresponding term in the 
expansion of the numerator. Further, we introduced the 
auxiliary function 



fi(x) = xip(l+x) 



(94) 



where ip(z) is the logarithmic derivative of the gamma 
function. The remaining terms contain 



K2 = gJ2&*[fi(tfc/u)] 



(95) 



and 



A, : , <( 2 [ 7 - HD/y)} + £ Re[f{(i&/i 



(96) 
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Whereas K\ and K 2 are independent of the cutoff, the 
last function diverges in the infinite bandwidth limit. For 
the average charge number we get likewise 



(n) = (n} c (l + (ffi) ) - (nK 1 + K 2 ) c . 



(97) 



Next, we determine {5nh) s where the time derivative 
of the charge number 



H > ^ = \ E (w«L°«<r A - H - 



kg a 



(98) 



arises from the tunneling Hamiltonian. Since h is of first 
order in the tunneling Hamiltonian and is off-diagonal 
in charge representation, we have to expand the density 
matrix up to first order, cf. Eq. (|90|), and one readily 
finds (Snn) s — 0(g 2 ). Now, the only term contributing 
in second order in the tunneling Hamiltonian reads 



{fin r (t)) s = — Re[tr{ 



e~ 0Ho he^ Hot ne~ 



£#oij 




-i(e+Sc)t\ t (99) 
' c 



where we have used that exp(iPL Pt)h = 
exp(j-H t)nexp(—j-H t) + 0(g). The Laplace transform 
can be readily performed leading to 



(*) = 



2(Sn 2 ) c 



?. sr 



<c + i C' nz 



C,C'=± 



o/3«C-»'C'fi2) 



^ z^ r m _ (trip, _ 



>„,!/£> 



^ [m- (iC'Zc - hz)/v][m - (*£'& + te)/i/] 



|5 c +tC'R«|/c 



(100) 



which depends on the electronic bandwidth D and shows 
a 1/z dependence for z 3> -D/Ti. The spectral density 
of the charge fluctuations now follows from Eqs. ( pa ) 
and (B9ft. Since the measurement device suppresses fre- 
quencies above a cutoff lu m <C .D/Ti, we may focus on 
|z| <C D/ft and find 



(*) 



.9 



2(6n 2 ) 



c 



C,C'=± 



ai 



-g c + i('hz 

- e P(t<~K'Kz) 



^Ea 



* K £c + C?i z ^ 



+2gz 



K=± 

ln(D/i/) - 7 
(<5n 2 ) c 



(101) 



which determines the noise spectrum for frequencies 
M < Z?/fi: 



5 n (w) = 2Re- 



C(0) 



(102) 



-jo; - 



-IOJ) 



While <j)(z) has been evaluated to first order in g the 
spectrum S n (w) contains terms of all orders in g. 



Two approximations can be considered: 
i) For large |w| 3> exp(—(8n 2 )c/g)D/h we may expand 
the denominator in Eq. (|102|) 



C(0) 



-iw + 4>(—iu>) 



(7(0) 



1 



</>(— iw) 



C(5 2 ) 



(103) 



and get the perturbative result which reads for |w| <C 

D/n 



S n (w) = 27rC / (0)<5(u;) 



^w 2 E ( ! 



C'/iw - £c 



cc=± *- / c 



in accordance with earlier findings 



0(g 2 ), (104) 

Here, the diverg- 
-iu) in Eq. (pip 



ing parts of C(0) in Eq. (96) and of i 
cancel leading to the finite contribution C-^(O) = C(0) 
(Ks)c independent of the cutoff I?. However, this ap- 
proximation is not valid at small frequencies and shows a 
1/lo 2 divergence that cannot describe the low frequency 
behavior correctly. In contrast, for large frequencies the 
resul t (104) is well behaved and merges with our result 

(Hi- 

ii) On the other hand, for small z -c nkBT/jh 
at large temperatures ksT ^> Ec and/or z <C 
ttEc exp(—PEc)/hln(/3Ec) at moderate to low temper- 
atures ksT~Ec we may replace <f)(z) by 0(0) + 4>'(0)z 
which is often referred to as Markovian approximation. 
Then 



C(z) 



C(0) 



z[l + 0'(O)]+0(O) 



(105) 



where <£'(0) = ~(K 3 ) c /{8n 2 ) C - In the limit 

2gln(D/u) <C (8n 2 )c we may expand 1 + (f>'(0) around 
g = and get 



S n (u) = Cf(0)- 



20(0) 



uJ- 



m 2 



(106) 



While the Markovian approximation fails to describe the 
high frequency behavior, at high temperatures the clas- 
sical result 



C m (uj) = 



k B TC 2G T /C 
e 2 uj 2 + (G T /C) 2 



(107) 



is recovered. He nce, at high temperatures the first order 
approximation (101) of the memory kernel <j)(t) becomes 
exact. 

Now, provided the experimental cutoff lu m is small 
enough for the Markovian approximation to be justified 
we have 

(Sn(t)Sn(0)) s M = i- f diuS n {u)F(u;/Lo M )e- wt (108) 
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with a cutoff function F(uo/lo m ) obeying F(Q) = 1. 
contour integration we get in the limit t — > 



Using 



(5n 2 ) M = Cf(0)F$(0)/w M ) = Cf(0) 



(109) 



Hence, for 0(0) <C oj m <C -D/fi. the measured variance is 
independent of the cutoff. 

To show the impact of finite tunneling conductance 
we compare in Fig. [f^ the variance fll09j ) of the average 
charge number in the Markovian approximation for a = 
(thin lines) and a = 1 (bold lines) for different temper- 
atures in dependence on the dimensionless gate voltage 
n g . Whereas tunneling amplifies the noise near n g = 0, 
in the vicinity of the degeneracy point the noise is sup- 
pressed by tunneling leading to the asymptotic behavior 



(6n 2 ) 



M 



g^n.{f3Ec/"K) + 1 + 7]. Since fluctuations 



of the charge are related to the linear conductance of the 
single electron transistor (SET) this behavior is in accor- 
dance with the observation that the linear conductance of 
the SET increases with the tunneling strength off the de- 
generacy point and on the other hand decreases directly 
at the degeneracy point, cf. Refs. |JLMB5|. 
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FIG. 15. Variance of the average charge number in depen- 
dence on the dimensionless gate voltage for different tempera- 
tures. Depicted are the results for a — (thin lines) compared 
with results for a — 1 (bold lines). 

Finally, we remark that directly at the degeneracy 
point 4>(z) becomes nonanalytic at zero temperatures for 
z — > and the Markovian approximation breaks down. 
Therefore, at lower temperatures fc^T <C Ec one has to 



IX. CONCLUSIONS 

In this article charge fluctuations of the single electron 
box were investigated by means of perturbation theory. 
Terms of third order in the phenomenological tunneling 
conductance were calculated analytically in the zero tem- 
perature limit. The predictions for the average charge 
number and the effective charging energy have been com- 
pared with Monte-Carlo data and renormalization group 
results. It has been shown that the perturbative treat- 
ment, in spite of its diverging behavior at the degeneracy 
point, leads to reliable results in a large part of the pa- 
rameter regime explored experimentally. An additional 
order in the perturbative series increases the range of va- 
lidity considerably, in particular for small gate voltages 
n g ss where PT covers the largest range of conduc- 
tance parameters a. When increasing the gate voltage 
the range shrinks continuously down to zero at the de- 
generacy point n g = |. In the vicinity of the degen- 
eracy point we have evaluated all graphs contributing to 
the leading logarithmic divergencies. A resummation was 
found to be insufficient to describe the behavior of the 
average charge number in the strong tunneling regime 
quantitatively. It turns out that nonleading logarithmic 
divergencies are essential, the summation of which re- 
mains an open problem. 

Further, the noise spectrum of the average charge num- 
ber has been investigated. In contrast to the average 
charge number this quantity depends on the electronic 
bandwidth D. Using a projection operator technique we 
have obtained an expression covering the perturbative 
as well as the semiclassical results. We find that in the 
Coulomb blockade region near n g = strong tunneling 
leads to an increase of the noise whereas at the degener- 
acy point n g = 3 the noise is suppressed. 
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use the full bandwidth dependent expression ( [10 0|) 



APPENDIX: THIRD ORDER RESULT 

In this appendix we give the explicit analytic result of the contribution in the third order in g to the ground state 
energy E . According to the transcendental functions appearing in the various terms we split the result into six 
contributions 



1G 



g 3 (u) = P3(w) + P2(w) + L3(u) + L2(u) + Ll(u) + i?(u) . 



(Al) 



Here P3 contains all terms with trilogarithms LI3 while PI lists those with dilogarithms Li2. Terms containing 
logarithms and no other transcendental hmctions are split into three types: In L3 and L2 expressions containing 
In and In are listed, respectively. Simple logarithms appear in LI and the remaining rational functions of u and 
constants are gathered in R. With the abbreviation 



these terms read 
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